B cells in the pneumococcus-infected lung are heterogeneous and require CD4+ T cell help including CD40L to become resident memory B cells.

Recovery from respiratory pneumococcal infections generates lung-localized protection against heterotypic bacteria, mediated by resident memory lymphocytes. Optimal protection in mice requires re-exposure to pneumococcus within days of initial infection. Serial surface marker phenotyping of B cell populations in a model of pneumococcal heterotypic immunity revealed that bacterial re-exposure stimulates the immediate accumulation of dynamic and heterogeneous populations of B cells in the lung, and is essential for the establishment of lung resident memory B (BRM) cells. The B cells in the early wave were activated, proliferating locally, and associated with both CD4+ T cells and CXCL13. Antagonist- and antibody-mediated interventions were implemented during this early timeframe to demonstrate that lymphocyte recirculation, CD4+ cells, and CD40 ligand (CD40L) signaling were all needed for lung BRM cell establishment, whereas CXCL13 signaling was not. While most prominent as aggregates in the loose connective tissue of bronchovascular bundles, morphometry and live lung imaging analyses showed that lung BRM cells were equally numerous as single cells dispersed throughout the alveolar septae. We propose that CD40L signaling from antigen-stimulated CD4+ T cells in the infected lung is critical to establishment of local BRM cells, which subsequently protect the airways and parenchyma against future potential infections.

Generation of human alveolar epithelial type I cells from pluripotent stem cells.

Alveolar epithelial type I cells (AT1s) line the gas exchange barrier of the distal lung and have been historically challenging to isolate or maintain in cell culture. Here, we engineer a human in vitro AT1 model system via directed differentiation of induced pluripotent stem cells (iPSCs). We use primary adult AT1 global transcriptomes to suggest benchmarks and pathways, such as Hippo-LATS-YAP/TAZ signaling, enriched in these cells. Next, we generate iPSC-derived alveolar epithelial type II cells (AT2s) and find that nuclear YAP signaling is sufficient to promote a broad transcriptomic shift from AT2 to AT1 gene programs. The resulting cells express a molecular, morphologic, and functional phenotype reminiscent of human AT1 cells, including the capacity to form a flat epithelial barrier producing characteristic extracellular matrix molecules and secreted ligands. Our results provide an in vitro model of human alveolar epithelial differentiation and a potential source of human AT1s.

Influenza Induces Lung Lymphangiogenesis Independent of YAP/TAZ Activity in Lymphatic Endothelial Cells.

The lymphatic system consists of a vessel network lined by specialized lymphatic endothelial cells (LECs) that are responsible for tissue fluid homeostasis and immune cell trafficking. The mechanisms for organ-specific LEC responses to environmental cues are not well understood. We found robust lymphangiogenesis during influenza A virus infection in the adult mouse lung. We show that the number of LECs increases 2-fold at 7 days post-influenza infection (dpi) and 3-fold at 21 dpi, and that lymphangiogenesis is preceded by lymphatic dilation. We also show that the expanded lymphatic network enhances fluid drainage to mediastinal lymph nodes. Using EdU labeling, we found that a significantly higher number of pulmonary LECs are proliferating at 7 dpi compared to LECs in homeostatic conditions. Lineage tracing during influenza indicates that new pulmonary LECs are derived from preexisting LECs rather than non-LEC progenitors. Lastly, using a conditional LEC-specific YAP/TAZ knockout model, we established that lymphangiogenesis, fluid transport and the immune response to influenza are independent of YAP/TAZ activity in LECs. These findings were unexpected, as they indicate that YAP/TAZ signaling is not crucial for these processes.

Total RNA sequencing reveals gene expression and microbial alterations shared by oral pre-malignant lesions and cancer.

Head and neck cancers are a complex malignancy comprising multiple anatomical sites, with cancer of the oral cavity ranking among the deadliest and the most disfiguring cancers globally. Oral cancer (OC) constitutes a subset of head and neck cancer cases, presenting primarily as tobacco- and alcohol-associated oral squamous cell carcinoma (OSCC), with a 5-year survival rate of ~ 65%, partly due to the lack of early detection and effective treatments. OSCC arises from premalignant lesions (PMLs) in the oral cavity through a multi-step series of clinical and histopathological stages, including varying degrees of epithelial dysplasia. To gain insights into the molecular mechanisms associated with the progression of PMLs to OSCC, we profiled the whole transcriptome of 66 human PMLs comprising leukoplakia with dysplasia and hyperkeratosis non-reactive (HkNR) pathologies, alongside healthy controls and OSCC. Our data revealed that PMLs were enriched in gene signatures associated with cellular plasticity, such as partial EMT (p-EMT) phenotypes, and with immune response. Integrated analyses of the host transcriptome and microbiome further highlighted a significant association between differential microbial abundance and PML pathway activity, suggesting a contribution of the oral microbiome toward PML evolution to OSCC. Collectively, this study reveals molecular processes associated with PML progression that may help early diagnosis and disease interception at an early stage.

Airway stem cell reconstitution by the transplantation of primary or pluripotent stem cell-derived basal cells.

Life-long reconstitution of a tissue's resident stem cell compartment with engrafted cells has the potential to durably replenish organ function. Here, we demonstrate the engraftment of the airway epithelial stem cell compartment via intra-airway transplantation of mouse or human primary and pluripotent stem cell (PSC)-derived airway basal cells (BCs). Murine primary or PSC-derived BCs transplanted into polidocanol-injured syngeneic recipients give rise for at least two years to progeny that stably display the morphologic, molecular, and functional phenotypes of airway epithelia. The engrafted basal-like cells retain extensive self-renewal potential, evident by the capacity to reconstitute the tracheal epithelium through seven generations of secondary transplantation. Using the same approach, human primary or PSC-derived BCs transplanted into NOD scid gamma (NSG) recipient mice similarly display multilineage airway epithelial differentiation in vivo. Our results may provide a step toward potential future syngeneic cell-based therapy for patients with diseases resulting from airway epithelial cell damage or dysfunction.

ß-catenin/CBP activation of mTORC1 signaling promotes partial epithelial-mesenchymal states in head and neck cancer.

Head and neck cancers, which include oral squamous cell carcinoma (OSCC) as a major subsite, exhibit cellular plasticity that includes features of an epithelial-mesenchymal transition (EMT), referred to as partial-EMT (p-EMT). To identify molecular mechanisms contributing to OSCC plasticity, we performed a multiphase analysis of single cell RNA sequencing (scRNAseq) data from human OSCC. This included a multiresolution characterization of cancer cell subgroups to identify pathways and cell states that are heterogeneously represented, followed by casual inference analysis to elucidate activating and inhibitory relationships between these pathways and cell states. This approach revealed signaling networks associated with hierarchical cell state transitions, which notably included an association between ß-catenin-driven CREB-binding protein (CBP) activity and mTORC1 signaling. This network was associated with subpopulations of cancer cells that were enriched for markers of the p-EMT state and poor patient survival. Functional analyses revealed that ß-catenin/CBP induced mTORC1 activity in part through the transcriptional regulation of a raptor-interacting protein, chaperonin containing TCP1 subunit 5 (CCT5). Inhibition of ß-catenin-CBP activity through the use of the orally active small molecule, E7386, reduced the expression of CCT5 and mTORC1 activity in vitro, and inhibited p-EMT-associated markers and tumor development in a murine model of OSCC. Our study highlights the use of multiresolution network analyses of scRNAseq data to identify targetable signals for therapeutic benefit, thus defining an underappreciated association between ß-catenin/CBP and mTORC1 signaling in head and neck cancer plasticity.

Structure learning for gene regulatory networks.

Inference of biological network structures is often performed on high-dimensional data, yet is hindered by the limited sample size of high throughput "omics" data typically available. To overcome this challenge, often referred to as the "small n, large p problem," we exploit known organizing principles of biological networks that are sparse, modular, and likely share a large portion of their underlying architecture. We present SHINE-Structure Learning for Hierarchical Networks-a framework for defining data-driven structural constraints and incorporating a shared learning paradigm for efficiently learning multiple Markov networks from high-dimensional data at large p/n ratios not previously feasible. We evaluated SHINE on Pan-Cancer data comprising 23 tumor types, and found that learned tumor-specific networks exhibit expected graph properties of real biological networks, recapture previously validated interactions, and recapitulate findings in literature. Application of SHINE to the analysis of subtype-specific breast cancer networks identified key genes and biological processes for tumor maintenance and survival as well as potential therapeutic targets for modulating known breast cancer disease genes.

Convergence of YAP/TAZ, TEAD and TP63 activity is associated with bronchial premalignant severity and progression.

BACKGROUND: Bronchial premalignant lesions (PMLs) are composed primarily of cells resembling basal epithelial cells of the airways, which through poorly understood mechanisms have the potential to progress to lung squamous cell carcinoma (LUSC). Despite ongoing efforts that have mapped gene expression and cell diversity across bronchial PML pathologies, signaling and transcriptional events driving malignancy are poorly understood. Evidence has suggested key roles for the Hippo pathway effectors YAP and TAZ and associated TEAD and TP63 transcription factor families in bronchial basal cell biology and LUSC. In this study we examine the functional association of YAP/TAZ, TEADs and TP63 in bronchial epithelial cells and PMLs. METHODS: We performed RNA-seq in primary human bronchial epithelial cells following small interfering RNA (siRNA)-mediated depletion of YAP/TAZ, TEADs or TP63, and combined these data with ChIP-seq analysis of these factors. Directly activated or repressed genes were identified and overlapping genes were profiled across gene expression data obtained from progressive or regressive human PMLs and across lung single cell RNA-seq data sets. RESULTS: Analysis of genes regulated by YAP/TAZ, TEADs, and TP63 in human bronchial epithelial cells revealed a converged transcriptional network that is strongly associated with the pathological progression of bronchial PMLs. Our observations suggest that YAP/TAZ-TEAD-TP63 associate to cooperatively promote basal epithelial cell proliferation and repress signals associated with interferon responses and immune cell communication. Directly repressed targets we identified include the MHC Class II transactivator CIITA, which is repressed in progressive PMLs and associates with adaptive immune responses in the lung. Our findings provide molecular insight into the control of gene expression events driving PML progression, including those contributing to immune evasion, offering potential new avenues for lung cancer interception. CONCLUSIONS: Our study identifies important gene regulatory functions for YAP/TAZ-TEAD-TP63 in the early stages of lung cancer development, which notably includes immune-suppressive roles, and suggest that an assessment of the activity of this transcriptional complex may offer a means to identify immune evasive bronchial PMLs and serve as a potential therapeutic target.

Polarity in respiratory development, homeostasis and disease.

The respiratory system is composed of a multitude of cells that organize to form complex branched airways that end in alveoli, which respectively function to guide air flow and mediate gas exchange with the bloodstream. The organization of the respiratory sytem relies on distinct forms of cell polarity, which guide lung morphogenesis and patterning in development and provide homeostatic barrier protection from microbes and toxins. The stability of lung alveoli, the luminal secretion of surfactants and mucus in the airways, and the coordinated motion of multiciliated cells that generate proximal fluid flow, are all critical functions regulated by cell polarity, with defects in polarity contributing to respiratory disease etiology. Here, we summarize the current knowledge of cell polarity in lung development and homeostasis, highlighting key roles for polarity in alveolar and airway epithelial function and outlining relationships with microbial infections and diseases, such as cancer.

Total RNA sequencing reveals gene expression and microbial alterations shared by oral pre-malignant lesions and cancer.

Head and neck cancers are a complex malignancy comprising multiple anatomical sites, with cancer of the oral cavity ranking among the deadliest and most disfiguring cancers globally. Oral cancer (OC) constitutes a subset of head and neck cancer cases, presenting primarily as tobacco-and alcohol-associated oral squamous cell carcinoma (OSCC), with a 5-year survival rate of ∻65%, partly due to the lack of early detection and effective treatments. OSCC arises from premalignant lesions (PMLs) in the oral cavity through a multi-step series of clinical and histopathological stages, including varying degrees of epithelial dysplasia. To gain insights into the molecular mechanisms associated with the progression of PMLs to OSCC, we profiled the whole transcriptome of 66 human PMLs comprising leukoplakia with dysplasia and hyperkeratosis non-reactive (HkNR) pathologies, alongside healthy controls and OSCC. Our data revealed that PMLs were enriched in gene signatures associated with cellular plasticity, such as partial EMT (p-EMT) phenotypes, and with immune response. Integrated analyses of the host transcriptome and microbiome further highlighted a significant association between differential microbial abundance and PML pathway activity, suggesting a contribution of the oral microbiome towards PML evolution to OSCC. Collectively, this study reveals molecular processes associated with PML progression that may help early diagnosis and disease interception at an early stage. AUTHOR SUMMARY: Patients harboring oral premalignant lesions (PMLs) have an increased risk of developing oral squamous cell carcinoma (OSCC), but the underlying mechanisms driving transformation of PMLs to OSCC remain poorly understood. In this study, Khan et al., analyzed a newly generated dataset of gene expression and microbial profiles of oral tissues from patients diagnosed with PMLs from differing histopathological groups, including hyperkeratosis not reactive ( HkNR ) and dysplasia, comparing these profiles with OSCC and normal oral mucosa. Significant similarities between PMLs and OSCC were observed, with PMLs manifesting several cancer hallmarks, including oncogenic and immune pathways. The study also demonstrates associations between the abundance of multiple microbial species and PML groups, suggesting a potential contribution of the oral microbiome to the early stages of OSCC development. The study offers insights into the nature of the molecular, cellular and microbial heterogeneity of oral PMLs and suggests that molecular and clinical refinement of PMLs may provide opportunities for early disease detection and interception.

Generation of human alveolar epithelial type I cells from pluripotent stem cells.

In the distal lung, alveolar epithelial type I cells (AT1s) comprise the vast majority of alveolar surface area and are uniquely flattened to allow the diffusion of oxygen into the capillaries. This structure along with a quiescent, terminally differentiated phenotype has made AT1s particularly challenging to isolate or maintain in cell culture. As a result, there is a lack of established models for the study of human AT1 biology, and in contrast to alveolar epithelial type II cells (AT2s), little is known about the mechanisms regulating their differentiation. Here we engineer a human in vitro AT1 model system through the directed differentiation of induced pluripotent stem cells (iPSC). We first define the global transcriptomes of primary adult human AT1s, suggesting gene-set benchmarks and pathways, such as Hippo-LATS-YAP/TAZ signaling, that are enriched in these cells. Next, we generate iPSC-derived AT2s (iAT2s) and find that activating nuclear YAP signaling is sufficient to promote a broad transcriptomic shift from AT2 to AT1 gene programs. The resulting cells express a molecular, morphologic, and functional phenotype reminiscent of human AT1 cells, including the capacity to form a flat epithelial barrier which produces characteristic extracellular matrix molecules and secreted ligands. Our results indicate a role for Hippo-LATS-YAP signaling in the differentiation of human AT1s and demonstrate the generation of viable AT1-like cells from iAT2s, providing an in vitro model of human alveolar epithelial differentiation and a potential source of human AT1s that until now have been challenging to viably obtain from patients.

Single Cell Transcriptomics of Fibrotic Lungs Unveils Aging-associated Alterations in Endothelial and Epithelial Cell Regeneration.

Lung regeneration deteriorates with aging leading to increased susceptibility to pathologic conditions, including fibrosis. Here, we investigated bleomycin-induced lung injury responses in young and aged mice at single-cell resolution to gain insights into the cellular and molecular contributions of aging to fibrosis. Analysis of 52,542 cells in young (8 weeks) and aged (72 weeks) mice identified 15 cellular clusters, many of which exhibited distinct injury responses that associated with age. We identified Pdgfra + alveolar fibroblasts as a major source of collagen expression following bleomycin challenge, with those from aged lungs exhibiting a more persistent activation compared to young ones. We also observed age-associated transcriptional abnormalities affecting lung progenitor cells, including ATII pneumocytes and general capillary (gCap) endothelial cells (ECs). Transcriptional analysis combined with lineage tracing identified a sub-population of gCap ECs marked by the expression of Tropomyosin Receptor Kinase B (TrkB) that appeared in bleomycin-injured lungs and accumulated with aging. This newly emerged TrkB + EC population expressed common gCap EC markers but also exhibited a distinct gene expression signature associated with aberrant YAP/TAZ signaling, mitochondrial dysfunction, and hypoxia. Finally, we defined ACKR1 + venous ECs that exclusively emerged in injured lungs of aged animals and were closely associated with areas of collagen deposition and inflammation. Immunostaining and FACS analysis of human IPF lungs demonstrated that ACKR1 + venous ECs were dominant cells within the fibrotic regions and accumulated in areas of myofibroblast aggregation. Together, these data provide high-resolution insights into the impact of aging on lung cell adaptability to injury responses.

Inactivation of LATS1/2 drives luminal-basal plasticity to initiate basal-like mammary carcinomas.

Basal-like breast cancers, an aggressive breast cancer subtype that has poor treatment options, are thought to arise from luminal mammary epithelial cells that undergo basal plasticity through poorly understood mechanisms. Using genetic mouse models and ex vivo primary organoid cultures, we show that conditional co-deletion of the LATS1 and LATS2 kinases, key effectors of Hippo pathway signaling, in mature mammary luminal epithelial cells promotes the development of Krt14 and Sox9-expressing basal-like carcinomas that metastasize over time. Genetic co-deletion experiments revealed that phenotypes resulting from the loss of LATS1/2 activity are dependent on the transcriptional regulators YAP/TAZ. Gene expression analyses of LATS1/2-deleted mammary epithelial cells notably revealed a transcriptional program that associates with human basal-like breast cancers. Our study demonstrates in vivo roles for the LATS1/2 kinases in mammary epithelial homeostasis and luminal-basal fate control and implicates signaling networks induced upon the loss of LATS1/2 activity in the development of basal-like breast cancer.

Virtual screening for small-molecule pathway regulators by image-profile matching.

Identifying the chemical regulators of biological pathways is a time-consuming bottleneck in developing therapeutics and research compounds. Typically, thousands to millions of candidate small molecules are tested in target-based biochemical screens or phenotypic cell-based screens, both expensive experiments customized to each disease. Here, our uncustomized, virtual, profile-based screening approach instead identifies compounds that match to pathways based on the phenotypic information in public cell image data, created using the Cell Painting assay. Our straightforward correlation-based computational strategy retrospectively uncovered the expected, known small-molecule regulators for 32% of positive-control gene queries. In prospective, discovery mode, we efficiently identified new compounds related to three query genes and validated them in subsequent gene-relevant assays, including compounds that phenocopy or pheno-oppose YAP1 overexpression and kill a Yap1-dependent sarcoma cell line. This image-profile-based approach could replace many customized labor- and resource-intensive screens and accelerate the discovery of biologically and therapeutically useful compounds.

Verteporfin protects against Th17 cell-mediated EAE independently of YAP inhibition.

The known YAP inhibitor verteporfin is capable of repressing IL-17A production in Th17 cells. However, this effect is mediated independently of YAP and can ameliorate Th17-mediated experimental autoimmune encephalomyelitis (EAE) upon in vivo administration. The data suggest verteprofin's mode of action for the design of novel therapeutic autoimmune disease intervention.

Basal extrusion of single-oncogenic mutant cells induces dome-like structures with altered microenvironments.

At the initial stage of carcinogenesis, oncogenic transformation occurs in single cells within epithelial layers. However, the behavior and fate of the newly emerging transformed cells remain enigmatic. Here, using originally established mouse models, we investigate the fate of RasV12-transformed cells that appear in a mosaic manner within epithelial tissues. In the lung bronchial epithelium, most majority of RasV12-transformed cells are apically extruded, whereas noneliminated RasV12 cells are often basally delaminated leading to various noncell-autonomous changes in surrounding environments; macrophages and activated fibroblasts are accumulated, and normal epithelial cells overlying RasV12 cells overproliferate and form a convex multilayer, which is termed a 'dome-like structure'. In addition, basally extruded RasV12 cells acquire certain features of epithelial-mesenchymal transition (EMT). Furthermore, the expression of COX-2 is profoundly elevated in RasV12 cells in dome-like structures, and treatment with the COX inhibitor ibuprofen suppresses the recruitment of activated fibroblasts and moderately diminishes the formation of dome-like structures. Therefore, basal extrusion of single-oncogenic mutant cells can induce a tumor microenvironment and EMT and generate characteristic precancerous lesions, providing molecular insights into the earlier steps of cancer development.

Inactivation of the Hippo tumor suppressor pathway promotes melanoma.

Melanoma is commonly driven by activating mutations in the MAP kinase BRAF; however, oncogenic BRAF alone is insufficient to promote melanomagenesis. Instead, its expression induces a transient proliferative burst that ultimately ceases with the development of benign nevi comprised of growth-arrested melanocytes. The tumor suppressive mechanisms that restrain nevus melanocyte proliferation remain poorly understood. Here we utilize cell and murine models to demonstrate that oncogenic BRAF leads to activation of the Hippo tumor suppressor pathway, both in melanocytes in vitro and nevus melanocytes in vivo. Mechanistically, we show that oncogenic BRAF promotes both ERK-dependent alterations in the actin cytoskeleton and whole-genome doubling events, which independently reduce RhoA activity to promote Hippo activation. We also demonstrate that functional impairment of the Hippo pathway enables oncogenic BRAF-expressing melanocytes to bypass nevus formation and rapidly form melanomas. Our data reveal that the Hippo pathway enforces the stable arrest of nevus melanocytes and represents a critical barrier to melanoma development.

Inhibition of LSD1 Attenuates Oral Cancer Development and Promotes Therapeutic Efficacy of Immune Checkpoint Blockade and YAP/TAZ Inhibition.

Lysine-specific demethylase 1 (LSD1) is a histone demethylase that contributes to the etiology of oral squamous cell carcinoma (OSCC) in part by promoting cancer stem cell phenotypes. The molecular signals regulated by LSD1, or acting with LSD1, are poorly understood, particularly in the development of OSSC. In this study, we show that conditional deletion of the Lsd1 gene or pharmacologic inhibition of LSD1 in the tongue epithelium leads to reduced development of OSCC following exposure to the tobacco carcinogen 4NQO. LSD1 inhibition attenuated proliferation and clonogenic survival and showed an additive effect when combined with the YAP inhibitor Verteporfin. Interestingly, LSD1 inhibition upregulated the expression of PD-L1, leading to immune checkpoint inhibitor therapy responses. IMPLICATIONS: Collectively, our studies reveal a critical role for LSD1 in OSCC development and identification of tumor growth targeting strategies that can be combined with LSD1 inhibition for improved therapeutic application.

Obesity-induced senescent macrophages activate a fibrotic transcriptional program in adipocyte progenitors.

Adipose tissue fibrosis is regulated by the chronic and progressive metabolic imbalance caused by differences in caloric intake and energy expenditure. By exploring the cellular heterogeneity within fibrotic adipose tissue, we demonstrate that early adipocyte progenitor cells expressing both platelet-derived growth factor receptor (PDGFR) α and ß are the major contributors to extracellular matrix deposition. We show that the fibrotic program is promoted by senescent macrophages. These macrophages were enriched in the fibrotic stroma and exhibit a distinct expression profile. Furthermore, we demonstrate that these cells display a blunted phagocytotic capacity and acquire a senescence-associated secretory phenotype. Finally, we determined that osteopontin, which was expressed by senescent macrophages in the fibrotic environment promoted progenitor cell proliferation, fibrotic gene expression, and inhibited adipogenesis. Our work reveals that obesity promotes macrophage senescence and provides a conceptual framework for the discovery of rational therapeutic targets for metabolic and inflammatory disease associated with obesity.

Epithelial LIF signaling limits apoptosis and lung injury during bacterial pneumonia.

During bacterial pneumonia, alveolar epithelial cells are critical for maintaining gas exchange and providing antimicrobial as well as pro-immune properties. We previously demonstrated that leukemia inhibitory factor (LIF), an IL-6 family cytokine, is produced by type II alveolar epithelial cells (ATII) and is critical for tissue protection during bacterial pneumonia. However, the target cells and mechanisms of LIF-mediated protection remain unknown. Here, we demonstrate that antibody-induced LIF blockade remodels the lung epithelial transcriptome in association with increased apoptosis. Based on these data, we performed pneumonia studies using a novel mouse model in which LIFR (the unique receptor for LIF) is absent in lung epithelium. Although LIFR is expressed on the surface of epithelial cells, its absence only minimally contributed to tissue protection during pneumonia. Single-cell RNA-sequencing (scRNAseq) was conducted to identify adult murine lung cell types most prominently expressing Lifr, revealing endothelial cells, mesenchymal cells, and ATIIs as major sources of Lifr. Sequencing data indicated that ATII cells were significantly impacted by pneumonia, with additional differences observed in response to LIF neutralization, including but not limited to gene programs related to cell death, injury, and inflammation. Overall, our data suggest that LIF signaling on epithelial cells alters responses in this cell type during pneumonia. However, our results also suggest separate and perhaps more prominent roles of LIFR in other cell types, such as endothelial cells or mesenchymal cells, which provide grounds for future investigation.